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This paper gives solutions for multilayer dielectrics that yield polariza- 
tion-independent operation at a given angle of incidence and at a given 
frequency. An exhaustive analysis is given for symmetric three-layer 
dielectrics surrounded by the same medium. Special solutions are also 
given for symmetric multilayer dielectrics with more than three layers 
surrounded by the same medium and for multilayer dielectrics surrounded 
by different media. Rules are presented for obtaining new polarization- 
independent solutions by cascading given solutions. Finally, for the purpose 
of demonstration and comparison of various solutions, different design 
examples are carried out for polarization-independent 3-dB beam splitters 
to operate at millimeter toavelengths. 

I. INTRODUCTION 

The state of polarization of a plane wave will, in general, be changed 
upon transmission through or reflection by a multilayer dielectric at 
oblique incidence. The reason for this depolarization is that the two 
eigenmodes of polarization, 2?-field perpendicular to the plane of inci- 
dence and i?-field parallel to the plane of incidence, generally have 
different reflection coefficients and different transmission coefficients. 
This is not desirable in many applications. For example, such depolari- 
zation could cause crosstalk in dual-polarization optical or quasi- 
optical components such as interference filters, 1 directional couplers 2-5 
and diplexers, 6 or mode conversion in beam-splitter-type hybrids 75 
in oversized waveguides. 

This paper provides solutions for multilayer dielectrics whose reflec- 
tion and transmission coefficients are independent of the state of 
polarization of the incident wave at a given angle of incidence and a 
given frequency. Compromising the degree of polarization indepen- 
dence to allow for a broader bandwidth of operation or for variations 
in the angle of incidence, even though desirable in many practical ap- 
plications, is not discussed. 
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The paper is primarily concerned with symmetric, multilayer di- 
electrics surrounded by the same medium. However, the case of differ- 
ent input and output media and antireflective coatings is also dis- 
cussed. Some of the results given in this paper have been previously 
reported by Kard, 8 Baumeister, 9 and Rabinovitch and Pagis. 10 Their 
work is referred to in appropriate places in the text. 

There are other types of structures, besides the ones reported in this 
paper, that can give polarization-independent operation at oblique 
incidence. For example, artificial or natural anisotropic dielectric 
layers or metallic-wire meshes with rectangular, rather than square, 
cells can, in principle, be constructed to achieve this effect. 

II. BASIC EQUATIONS 

Consider n adjacent isotropic dielectric layers of infinite transverse 
extent having relative dielectric constants of e h a, ■•-,«« and uniform 
thicknesses of di,d 2 , ■ ■ • , d n inserted in a medium of a relative dielectric 
constant e r . Define 



- e,/« 



i = 1, 2, 



(1) 



as the normalized dielectric constants of the layers. Let a uniform 
plane wave propagating in the medium be incident at an angle 6 on the 
layers as shown in Fig. 1. The incident, transmitted and reflected waves 
will have the same state of polarization if the i£-field is perpendicular 




Fig. 1 — n-layer dielectric inserted in a medium of a relative dielectric constant « r . 
The subscript h refers to the /t-mode (i?-field perpendicular to plane of incidence), 
and e refers to the e-mode (JS'-field parallel to plane of incidence). 
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to the plane of incidence (h-mode), or if the E-field is parallel to the 
plane of incidence (e-mode). To calculate the transmission and reflec- 
tion coefficients for the h- and e-modes, one can represent each layer 
by an equivalent transmission line. 11,12 The impedance of the line 
representing the ith layer normalized to that of the line representing 
the medium in which the layers are inserted is 



Zi = Z ih = cos d/im - sin 2 0)* 
for the /i-mode, and 

Zi = Z ie = (m - sin 2 0) */(«.- cos 0) 
for the e-mode. The corresponding electrical length is 

<t> { = {2ir/\ )d i 47r{K i - sin 2 0)i 



(2a) 



(3a) 



(4a) 



for both modes; X being the free-space wavelength. If k, < sin 2 0, i.e., 
the total reflection condition is satisfied at the boundary of the ith 
layer, then the square root of (k, — sin 2 0) is imaginary and one can 
write 

Z ih = jX ih , (2b) 



Zie — jXie, 

0t = — j<Xi, 



(3b) 

(4b) 



where X«, Xu, and a, are positive, real quantities. In this case, cou- 
pling is accomplished across that layer by an evanescent wave. 
The normalized ABCD matrix corresponding to the ith. layer is 



M,= 



COS tf>i 



j sin 4>iZ t 



j sin 4n/Zi cos 0; 

The overall, normalized ABCD matrix of the n layers is 
A jB' 



JC 



D 



= M 1 M 2 .--M„. 



(5) 



(6) 



The reason for the j in the off-diagonal terms of the matrix of eq. (6) 
is that A, B, C, and D will be real quantities if the system is lossless. 
Henceforth, this will be assumed to be the case. 
It is worth noting that reciprocity requires that 

AD + BC = 1, (7) 

and that for a symmetric structure 

A = D. (8) 
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The overall field transmission coefficient, t, and field reflection co- 
efficient, r, are given by 

t m Et/Ei = 2/lA + D + j(B + C)], (9) 

r = &/JB, = [A - D + j(B - C)]/[A + Z) + j(5 + C)]. (10) 

A subscript of h or e should be added to all the variables in eqs. (5) 
through (10) to denote the particular mode being considered. 

III. CLASSIFICATION OF THE SOLUTIONS 

Our goal is to find solutions for the k's, the d's, and to achieve 
polarization-independent operation. In applications where one is 
restricted to independent h- and e-modes, or to unpolarized waves, it 
is sufficient that 

\h\ = \te\, \r h \ = \r e \. (11) 

To satisfy these conditions, eqs. (7), (9), and (10) require that 

Al + B\ + Cl + Dl = A\ + B\ + CI + Dl (12) 

For symmetric structures, eqs. (7) and (8) reduce eq. (12) to 

B h -C h = ±(B e -C e ). (13) 

To preserve the state of polarization of an arbitrarily polarized 
incident wave, condition (11) is not sufficient. Instead, it is necessary 
that 

t k = ±t e , r h = ±r e . (14) 

The sense of polarization (i.e., the right- or left-handedness) of the 
transmitted wave will be unchanged if fa = -\-t e and will be reversed 
if h = —t e . On the other hand, the sense of polarization of the reflected 
wave will be unchanged if tk = —r e and will be reversed if ru = -\-r e . 
The latter case is similar to reflection by a perfect, conducting plane. 
The four combinations of signs in eq. (11) are listed in Table I where 
they are referred to as Cases 1 through 4. The required relations 

Table I — Four cases for polarization-independent operation 



\Case 


1 


2 


3 


4 


h = 


t. 


t. 


-u 


-u 


r* = 


r. 


— r« 


r. 


—r. 


A h = 


A. 


D e 


-A e 


-D. 


B h = 


B. 


c. 


-B e 


-c. 


C h = 


c. 


B. 


-c. 


-B e 


D h = 


D e 


A e 


-D. 


-A. 
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between A h , B h , C h , D h and A e , B e , C e , D e to satisfy each of these cases 
can be deduced from eqs. (9) and (10) and are also given in Table I. 
In addition to classifying the solutions in the manner of Table I, 
another classification based on the numerical values of the k's is useful. 
To be concise, let K mia be the smallest value of the k's and let it be at 
the rath layer, i.e., 



Kmin — Him (Ki, Ki 



j K n ) Km- 



(15) 



If K m j„ ^ 1, then none of the k's is smaller than unity and the solution 
can be realized by dielectric layers, e» ^ 1, inserted in air, e r = 1. This 
will be referred to as a Type A solution. On the other hand, if K min < 1, 
the surrounding medium can no longer be air since this would require 
that at least e m be less than unity (recall that m is the layer with the 
smallest value of k) ; thus, a dielectric-prism realization as shown in 
Fig. 2 is necessary. In this case, the wth layer can be an air separation 
e m = 1, and thus, from eq. (1) 



€, = Kj/Kmia = 1> e r = lAmin > 1- 



(16) 



The prism realization will be referred to as a Type B solution when 
1 > K min ^ sin 2 0, and a Type C solution when /c min < sin 2 0. The dis- 




MATCHING LAYERS 



Fig. 2 — A prism realization having the mth layer (the layer with the smallest 
dielectric constant) made of an air separation. Note that k„ = Kmia < 1- Matching 
layers are required at the external surfaces of the prisms. 
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Table II — The three types of solutions 



Type 


A 


B 


C 


Condition 


itmin ^ 1 


Sin 2 ^ K m in < 1 


« m in < sin 2 B 


Description 


Dielectric slabs in 
air (< r = 1 ) is 
possible. 


Prism realization 

(«r > 1) is 

necessary. No 

evanescent 

waves. 


Prism realization 
(e r > 1) is neces- 
sary. Evanescent 
waves exist. 



tinction is made because in a Type C solution, the condition for total 
reflection is satisfied at the boundary of the wth layer and coupling 
across that layer is accomplished by an evanescent wave, as mentioned 
in Section II. This is not the case for a Type B solution. Table II 
summarizes the conditions for the three types of solution. 

Dielectric-prism devices have been used for many optical 13 ' 1 * and 
quasi-optical 2-5 - 15-18 applications. When the total reflection condition is 
satisfied, Type C realization, the device is often said to be of the frus- 
trated-total-reflection type. 

Fresnel reflections at the outer surfaces of the prisms might cause 
undesirable effects. 19 However, this can be avoided by matching each 
of these surfaces by the familiar X/4 transformer, or by other 
means. 2 ' 20-22 In any case, for the purpose of this paper, the method of 
matching should be independent of the state of polarization of the 
incident wave. Thus, Brewster angle matching 1718 would not be suit- 
able here. 

In the remainder of the paper, a combined reference to the four 
cases of Table I and the three types of Table II is adopted. For ex- 
ample, Case 1A refers to Case 1 of Table I and Type A of Table II. 

IV. SINGLE-LAYER SOLUTIONS 

It can be shown by using the procedure outlined in Section II that 
no solution exists for a single layer to satisfy any of the four cases 
indicated in Table I. However, it is possible to satisfy the weaker con- 
dition given in eq. (11) by requiring that 

Kl = ei / Cr = sin 2 0/(1 + cos 2 6), (17) 

where the symbols are defined in Fig. 1 with n = 1. Substituting eq. 
(12) into eqs. (2) through (10), one obtains 

k = C, r h = r*, (18) 

where the asterisk denotes complex conjugation. 
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It is clear from eq. (17) that «i < sin 2 and, hence, from Table II, 
the solution is of Type C, and a prism realization is necessary. 

It is interesting to note that condition (17) is independent of the 
thickness, d\, of the layer. Thus, varying d\ will change the reflection 
and transmission coefficients without affecting the polarization inde- 
pendence. If R = |r| 2 is the required power reflection coefficient, then 
eqs. (17) and (2) through (10), with n = 1, give 

d x = [X /(2»VT 1 )][(1 + *i)/(l - n)]* 

X sinh- 1 (20^/(1 - B)]*/(l + «0I- (19) 

A continuous variation of d x , and hence of R, can be accomplished 
most conveniently if the layer is simply an air separation, e x = 1. In 
this case, if the angle of incidence is = 45°, which is convenient for 
mechanical considerations, eq. (17) requires the dielectric constant 
of the prisms to be e r = 3. On the other hand, if the prisms were made 
of polystyrene, e r = 2.54, eq. (17) requires that = 48.73°. 

V. SYMMETRIC THREE-LAYER SOLUTIONS 

It was mentioned in the previous section that no single-layer solution 
exists that satisfies any of the four cases of Table I. The same is also 
true for two layers. On the other hand, one can satisfy any of the four 
cases with three or more layers. In this section, we give all the solutions 
for symmetric three-layer structures, i.e., 

k 3 = k 1} d 3 = rfi. (20) 

Even though unsymmetric three-layer structures can have polariza- 
tion-independent solutions, they will not be discussed here because of 
the algebraic complexities involved. 

Using the terminology of Section II, it can be shown that the 
normalized A BCD parameters of a symmetric three-layer structure 
are given by 

A = D = cos 2 cos (200 - h[Zi/Zi + Z«JZ{] sin <f> 2 sin (200 (21) 

B = Z\ cos 02 sin (200 + Z» sin 2 cos 2 0i 

- [Z 2 /Z,] sin 2 sin 2 0i (22) 

C = Zf 1 cos 02 sin (200 + Z7 l sin 2 cos 2 0i 

- [Z./Z?] sin 2 sin 2 0l (23) 

Combining these equations with eqs. (2) through (4), (9), and (10), 
one can obtain solutions for the cases listed in Tables I and II. In 
Cases 3 and 4 of Table I, this process involves considerable algebraic 
manipulation. The results are given below. For simplicity, expressions 
for 0i (or txi) and 0j (or <x>) are given instead of expressions for d x and 
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d>. Also, we define 

Pi = (k, - sin 2 0)S, i = 1,2 (24a) 

g f = (sin 2 - *,)*, i= 1,2 (24b) 

/(«, y) = 2xy - (x 4- y) sin 2 5, (25) 

# = |r»|» = |r e | 2 , (26) 

S m [fl/(l - R)y, (27) 

U = [2Vfi/(l + Vfl)]», (28) 

V = (1 - Vfi)/(1 + Vfl). (29) 

In each case, necessary upper and/or lower bounds on k\, k 2 , 0, and R 
are given. Because of the algebraic complexity involved, the bounds on 

Ki and K 2 in Cases 3 and 4 are looser than those given in Cases 1 and 2. 

Case 1 A : 

0! = (2m + l)x/2, m = 0,1,2, •■-. (30) 

K , = K f/[l - ( Kl - l)Vcos 2 62. (31) 

sin 2 = ±2S/[V7»Ai - «i/V7 2 ]. (32) 

1 + [/ cos 6 g ki < 1 + cos 6. (33) 

(1 + U cos 0) 2 /F g k 2 < «. (34) 

This solution exists for all R, 6. 

Case IB: 

0i, K-i, </»2 are the same as in eqs. (30) through (32). 

sin 2 6 < ki ^ 1 - U cos 6. (35) 

sin 2 e < ko g> (i - u cos ey/v. (36) 

1 < KiA 2 ^ max [2/(1 + sin 6), V/(l - U cos 0)2. (37) 
This solution exists if and only if 

sin 6 < V\ (38a) 

or, equivalently, 

R < [cos 2 0/(1 + sin 2 0)] 2 . (38b) 

Case 1C: 

No solution exists. 
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Case 2 A : 

0x = (2m + l)ir/2, m = 0, 1, 2, • ••. (39) 

K2 = kt. (40) 

sin a = ±25 cos 0pfp 2 /[(*i - 1) sin 0] 2 . (41) 

*!>«■ (>1), (42) 

k 2 >*2 (>1), (43) 

where n a is the root larger than unity of the quartic equation 

(W - l)* 4 + 2(2 - W sin 2 0)k 3 - (6 + W sin 2 cos 2 0)* 2 

+ 2(2 + W sin 4 0)k - (1 + W sin 6 0) = 0, (44) 

where 

W = [25 cos 0/sin 2 0] 2 . (45) 

This solution exists if and only if 

sin e > U, (46a) 

or, equivalently, 

R < [sin 2 0/(1 + cos 2 0)] 2 . (46b) 

Case 2B : 

0i, k 2 , 2 are the same as in eqs. (39) through (41). 

sin0 < ki ^ K b (<1), (47) 

sin 2 0<K, ^ kI (<1), (48) 

where Kb is the positive root less than unity of the quartic equation (44). 
This solution exists for all R, 0. 

Case 2C: 

0i, k 2 are the same as in eqs. (39) and (40). 

sinh a 2 = 2S cos 0p?g 2 /[(l - *i) sin 0] 2 . (49) 

sin 2 < ki < sin 0. (50) 

sin 4 < k 2 < sin 2 0. (51) 

This solution exists for all R, 0. It has been previously found by 
Kard. 8 

Cases 3 A, SB: 

No solutions exist. 
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Case 3C : 

cos (200 = f(K h 1)/(ki, <c 2 )/[(k! - 1)(k-i - k 2 ) sin 4 5]. (52) 

• l, _ 2V^g 2 /(Kl, 1) ,-,v 

sinna 2 - (1 + Ki)(ki _ Ka) sin2 fl |- r _ (K2 + J2) C0S 2 fjjv W 

S = Vm(ic 2 - X) cos 0/17 - fo + X 2 ) cos 2 5]», (54) 

where 

X^ /0a,/c 2 )/[(l +K 1 ) S in 2 ^] ) (55) 

Y m k 2 (1 - K2 )(ki - 1)«/(1 + Kl). (56) 

ki > (1+3 cos 2 6 + [8 cos 2 5(1 + cos 2 0)]»}/sin 2 5. (57) 

sin 2 5 cos 2 5/(1 + cos 2 5) ^ k 2 ^ sin 2 5/(1 + cos 2 5). (58) 

This solution exists if and only if 

sin 5 > (2/S)*{[(l + 32/27S 2 )* + 1]* 

- [(1 + 32/27S 2 )* - 1]*}*, (59a) 

or, equivalently, 

R > 8 cos 2 5/(8 cos 2 5 + sin 6 5). (59b) 



Cases 4 A, 4B: 




No solutions exist. 




Case 4C: 




cosh2a 1 = (1 + k0(k 2 + Kj)/[(1 - ki)(«« - «i)]. 


(60) 


2Vk~i* 2 (1 + Ki)gip 2 


(61) 


tal102 /(Kl,K 2 )[(l+K 2 )W + K 2 )]i 


5 = V^(k 2 - X)/[cos 5 9i (k 2 + X 2 )*], 


(62) 


where X is given in eq. (55). 




«i < sin 2 5. 


(63) 


*c 2 > sin 2 5. 


(64) 


This solution exists for all R, 5. 




Discussion : 





The above exhausts all possible solutions for a symmetric three-layer 
dielectric with no polarization dependence. The choice of the par- 
ticular solution to implement in a given problem depends on the desired 
values of R and 5, the available dielectric materials, and whether 
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having dielectric layers surrounded by air (Type A solution) or a 
prism realization (Type B or C solution) is more suitable. 

Cases 1A and 2A, which are the only possible Type A solutions, 
require large values for the dielectric constants and/or the angle of 
incidence to achieve moderate or large values of R. This point is 
illustrated for R = 0.5 in the design examples given in Section VIII. 

In Case 2, the required relation between *i and k 2 given in eq. (40) 
has an advantage for the prism realizations in Cases 2B and 2C. In 
these two cases, since «•> < «i < 1, eqs. (16) and (40) give e 2 = 1, 
ei = ki/k 2 , and e r = 1/kj = *?■ Thus, ei can also be used as a quarter- 
wave transformer to match the outer surfaces of the prisms. 

As in the single-layer prism solution given in Section IV, one can 
vary R in Cases 1 and 2 without affecting the polarization independence 
by simply changing d>. This is particularly convenient in the prism 
realizations in Cases IB, 2B, and 2C with e 2 = 1 (air) since, in this 
case, d> can be changed continuously. Case 2C is unique in that one 
can obtain any value of R from zero to unity by varying d-> from zero 
to infinity. In Cases 1 A, IB, 2A, and 2B, varying d 2 results in a periodic 
variation of R between zero and a maximum value less than unity. 
Varying rii or d> in Cases 3 or 4 results in a polarization-dependent 
operation. 

Cases 1A, 2B, 2C, and 4C have solutions for all R and 0. This is not 
true for the remaining three cases, IB, 2A, and 3C, as indicated by eqs. 
(38), (46), and (59). However, one should note that the numerical 
values of e r , ei, and e 2 , and not necessarily the limitations on R and 6, 
often determine whether or not any particular case can be used in 
practice. 

In Case 4C, eq. (46) indicates that one can have tc 2 = 1. In this case, 
given that k x < 1 from eq. (63), eq. (16) gives d = 1 and c r = « 2 
= 1/jci. Thus, the prism solution can be realized by using one dielectric 
material with air separations. This is not true for any of the other cases. 

VI. SYMMETRIC MULTILAYER SOLUTIONS 

When there are more than three layers, an exhaustive analysis 
similar to that given in the previous section becomes quite involved. 
However, two relatively simple classes of solutions for a symmetric 
multilayer dielectric with no polarization dependence can be obtained 
by generalizing Cases 1 and 2 of the symmetric three-layer dielectric. 
First, one should note that symmetry requires that the number of 
layers n be odd, i.e., 

n = 2l + l, 1 = 1,2, .... (65) 
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Thus, from symmetry 

Ki = Kn+t-i, di = dn+l-i, i — 1, 2, • • • , I. (66) 

Let us assume that the electrical lengths of the first, and thus also 
the last, I layers are odd multiples of 7r/2; i.e., 

t - = (2m, + l)x/2; m, = 0, 1, 2, ■ ■ - ; »- 1,2, •••,/. (67) 

While this is a necessary condition to obtain solutions for Cases 1 and 2 
of a symmetric three-layer dielectric, as indicated in the previous 
section, eqs. (30) and (39), it is not necessary in general. However, 
this assumption is employed here because it simplifies the analysis 
considerably. 

Define the function (ji of the vector x = [x\, x 2 , • ■ -, xi + i} as 

j7,(x) - xiiv )i \ n *?)/( n *?). (68) 

Thus, g\(x) = x"\/x->, <72 (x) = x\xz/xl, • ■ • , etc. Each x, can assume the 
value of Zih, Zu, Ki, or p, (which will be defined later). From eqs. (66) 
and (67) and Section II, it can be shown that the normalized ABCD 
parameters of the symmetric multilayer dielectric are given by 

A = D = (-l)'cos(0, + i), (69) 

B = (-l)<sin(0, +1 )0,(Z), (70) 

C = (-l)'sin(0 /+1 )/^(Z). (71) 

Since <f>i + i is the same for both modes, eq. (69) shows that Ah = Dh 
= A e = D e - Thus, from Table I, the above equations can give solutions 
for Cases 1 and 2 provided that the conditions for the B's and the C's 
are met. Before satisfying these conditions, however, we note from eqs. 
(2) and (3) that 

ZaZu = l/Ki (72) 

and define 

Vi = Zi e /Z ih = ( Ki - sin 2 B)/(Ki cos 2 6). (73) 

As Ki assumes the increasing values of zero, sin 2 6, 1, and + o°, then vt 
assumes the increasing values of — °o, 0, 1, and sec 2 6, respectively. 

From eqs. (9), (10), and (65) through (73), and from Table I, one 
can obtain solutions for Cases 1 (B h = B e , Ch = C e ), and 2 (B h = C, 
Ch = B e )- The results are given below. Because of algebraic complexity, 
each case has not been subdivided into Type A, B, and C solutions. 
Also, for the same reason, no bounds are given on the /c's, 6, or R. 

Case 1 : 

gi(v) = 1. (74) 

sin<^ +1 = ±2S/[M*)y - Ode)]"*}. (75) 
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Case 2 : 

gM = 1. (76) 

sin 0/ + , = ±2S/|[r//(v)]5 - [ f y,(v)]-M. (77) 

In eqs. (75) and (77), S is the same function of R defined in eq. (27). 
Note that the roles of the k's and the v's are interchanged in the two 
cases. 

For a solution of any particular case to be valid, the value of sin </>j + i 
should either be real with a magnitude not exceeding unity, which 
leads to a Type A or a Type B solution, or be imaginary, which leads 
to a Type C solution. The latter solution is not possible in Case 1 be- 
cause the k's are real positive quantities, and thus the right-hand side 
of eq. (75) is always real. On the other hand, since the v's, and in par- 
ticular vi+i, can be negative, a Type C solution is possible in Case 2. 

As in the previous section, solutions for Cases 1A, 2B, and 2C in the 
present section are possible for all values of R and 0. 

In Cases 2B and 2C, when ki+i does not exceed k\, k 2 , • • • , ki, or unity, 
eq. (16) shows that one can have a prism realization with ej+i = 1, 
e t = ki/ki+i, ■ • •, «i ~ ki/ki+i, and e r = l/«i+i. Using eq. (76), one can 
show that / quarter-wave layers with dielectric constants ei, e 2 , • • • , u 
will match the outer surfaces of the prisms. 

In practice, it is desirable to construct the multilayer dielectric 
reflector with the least number of different dielectrics. Thus, let us 
assume that two dielectrics with normalized dielectric constants k x and 
k-2 are used alternatively; i.e., 

Ki = ki, i odd, (78a) 

Ki = k-2, i even. (78b) 

Combining this assumption with those in eqs. (65) through (67), the 
general solutions given in eqs. (74) through (77) simplify considerably. 
Because of its practical importance, we restrict ourselves to Type A 
solution, where k\ and k 2 are larger than unity. No solution exists if 
either ki or k-< is equal to unity. 

Case 1A : 

k-2 = sin 2 0/ 1 1 - [(/ci - sin 2 0)Ai]< l+1 >" (cos*)" 2 "}- (79) 

sin /tI = ±2S/[(k 2 A<+ 1 )« - (id +l /4)*]. (80) 

K , < sin 2 0/[l - (cos0) 2 ' ( ' +1) ]. (81) 



,'+i /J 



4 < V. (82) 



This solution exists for all R, 6. 
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Case 2 A : 

Ki = K ? +1)/l (83) 

on**, = ±2S/ZW l /W - MM")*]. (84) 

*,!,/,/+ ' < 7. (85) 

This solution exists if and only if 

sin 6 > U, (86a) 

or, equivalently, 

R < [sin 2 9/(1 + cos 2 6) J. (86b) 

The quantities S, U, and F used in the above equations are defined 
in eqs. (27) through (29). Of course, when I = 1, i.e., n = 3, the above 
solutions reduce to Cases 1A and 2A given in the previous section for 
symmetric three-layer dielectric. 

VII. CASCADING SOLUTIONS 

It may happen that some required values of the reflection coefficient 
cannot be obtained by any of the solutions given in the previous two 
sections because of limitations on the values of the dielectric constant 
available in practice, or because the resulting angle of incidence is not 
suitable for a particular application. In this case, the desired multi- 
layer dielectric may be obtained by cascading two or more solutions, 
provided that they have the same a r and the same 8. 

We now proceed to find the rules for proper cascading to maintain 
polarization-independent operation. Let the two solutions to be cas- 
caded be labeled x and y. They are not necessarily symmetric or 
identical. Let = 2ird cos 0/X be the electric length between x and y 
(see Fig. 3). Dropping the subscripts h and e, which denote the par- 
ticular mode being considered, we define r x , r„ t x , r y , r y , and t y as the 
reflection and transmission coefficients of x and y in the directions 
specified in the figure. The transmission coefficients t x and t v do not 
depend on the direction of propagation because of reciprocity and 
because the medium is the same on both sides of each of x and y. 
Because of losslessness, the difference between r« and r' x or between 
r„ and r'„ is a phase factor. Of course, if x is symmetric then r x = r x , and 
if y is symmetric then ?•„ = r' y . 

Using the method of successive reflections, 23 it can be shown that for 
a wave incident on the x side of the cascade, the overall field trans- 
mission coefficient, t, and field reflection coefficient, r, are given by 

t = t x t y exp (-J0)/[1 - r x r y exp (- j2<i>)~\, (87) 

v = r' x + tlr y exp (- J20)/[l - r x r y exp (-./20)], (88) 
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Fig. 3— Cascading of two polarization-independent solutions x and y to obtain an 
overall polarization-independent solution. 

where we recall that a subscript of h or e has been dropped from all the 
t's and r's. 

Since x and y are polarization-independent solutions, then, from eq. 
(14) or Table I, r xh r yh = ±r xe r ye . Only the positive sign is acceptable, 
otherwise the magnitude of the denominators in eqs. (87) and (88) 
would be different for the h- and e-modes, which will result in a polari- 
zation-dependent operation for the cascade. Thus, the sign of r e /r h 
should be the same for all the cascaded sections; i.e., one can only have 
a combination of Cases 1 and 3 or Cases 2 and 4. 

Let the abbreviation i*j = k mean: "Cascading Cases i and j gives 
Case k." It can be shown from eqs. (87) and (88) and Table I that 

1*1 = 3*3 = l, 1*3 = 3*1 = 3, (89a) 

2*2 = 4*4 = 2, 2*4 = 4*2 = 4. (89b) 

Let the required overall power reflection coefficient of the cascade be 

R = |r| 2 , (90) 

and let it be realized by two identical cascaded sections each having a 
power reflection coefficient R x . Then, one can write 

r x = r u = ft***., (91) 
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where <j> x is the phase angle of r x . This can be found given the e's, d's, 
and by using eqs. (1) through (10). Substituting eq. (91) in eqs. (87) 
and (88), and denning 

* = 4>* ~ 4>, (92) 

one obtains 

B = [I + (1 - R x )*/(4R x am*f)Jr l . (93) 

Thus, for a given value of R x , one can obtain any value of R satisfying 

S R ^ 4/2,/ (1 + R x ) 2 . (94) 

Note that the upper bound on R is larger than R x . Conversely, for a 
given value of R, one can use any value of R x satisfying 

[1 - (1 - R^J/R £R X £1. (95) 

The higher the value used for R x , the smaller is the overall bandwidth. 

VIII. DIFFERENT INPUT AND OUTPUT MEDIA, AND ANTIREFLECTION 
COATINGS 

Often, especially at optical wavelengths, the multilayer dielectric is 
supported by a substrate having a relative dielectric constant e„ which 
is, in general, different from that of the input medium, c r . In this case, 
general analysis of polarization-independent solutions is complicated. 
However, it can be shown that special solutions similar to those for the 
symmetric (21 + l)-layer dielectric given in Section VI can be easily 
obtained. For this purpose, let 6, as before, be the angle of incidence 
in the e r medium. Assume that there are J( = l) layers with relative 
dielectric constants ei, e->, • • • , a between the e r and e s media. The 
normalized dielectric constants «ci, k 2 , • • •, Kt are defined by eq. (1). 
In addition, we define 

ki+i = e s /e r - (96) 

Further, as was done in Section VI, assume that each of the I layers is 
effectively a quarter-wave thick at the given angle of incidence, i.e., 

0, = (2m,- + 1)t/2; rm = 0, 1, 2, • • -, i= 1,2, ■■ ■, I. (97) 

With the above definitions and assumptions, it can be shown from 
eqs. (1) through (6) that the normalized input impedance as seen from 
the e r medium is given by 

Z in = gi (Z), (98) 

where Z = \Z\, Z», ■ • • , Zi +i ] and the function gi is defined in eq. (68). 
Thus, the reflection coefficient in the e r medium is given by 

r = (Z in - l)/(Z in + 1). (99) 
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A subscript of h or e, depending on the particular mode being con- 
sidered, has been dropped from all the variables in eqs. (98) and (99). 
It is clear that if kj+i = eje r > sin 2 0, then Z i+h Z in and r are all 
real quantities. In this case, it can be shown that the reflection coeffi- 
cient r, in the e s medium, and the transmission coefficient / across the 
I layers, are given by 

r,. = (-1)'+V, (100) 

and 

t = (_/)<(_!) 2 »«(1 _ r 2 )>, (101) 

where the m,-'s are the integers defined in eq. (97). Since the impedances 
of the input and output media are different, t is denned here as the 
scattering transmission coefficient, which is usually referred to as S v >. 
Using the same notation, r and r„ could be referred to as Su and S->», 
respectively. 

The coefficients r s and t can not be defined if the total reflection con- 
dition, ki+i s e„/«r < sin 2 0, is satisfied. In this case, Zj+i is imaginary 
and a nonpropagating evanescent wave exists in the e, medium. It 
also follows that Z in is imaginary, and hence, from eq. (99), \r\ = 1. 

It can be seen from eq. (99) that for r e = r h (Case 1), Z in , e = Z*»,»; 
and that for r e = -r h (Case 2), Z in , e = \/Z in , h . Thus, combining eqs. 
(98), (68), (72), and (73), and defining R = \r e \ 2 = \ r h | 2 as the power 
reflection coefficient, one obtains the following solutions. 

Case 1 : 

gi(v) = 1. (102) 

R = {[>(*)]* - U7{|>(k)]» + IP- ( 103 > 

c,,( K ) = 1. (104) 

R= {{.(jMJ- i}V{[^(v)]* + i} 2 . d05) 

The above solution of Case 2 has been previously obtained by 
Baumeister 9 and by Rabinovitch and Pagis. 10 Both papers, however, 
overlooked the solution of Case 1. 

It is observed that eqs. (102) and (104) are identical to eqs. (74) and 
(76) for the symmetric (2/ + l)-layer dielectric. This should not be 
surprising since the latter can be obtained by cascading, back to back, 
two identical solutions of the type described in the present section. 
In fact, this is the reason that the conditions for polarization-inde- 
pendent operation given in Sections V and VI for Cases 1 and 2 are 
independent of the thickness of the middle layer. 
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Case 2 : 



It is clear from eqs. (102) through (105) that if 

9i(*) = 9i(y) = 1, (106) 

one obtains R = 0. Thus, eq. (106) gives the conditions for polariza- 
tion-independent antireflection coatings at inclined incidence. It has 
been previously reported by Bauneister. 9 Rabinovitch and Pagis 10 
have incorrectly stated that no such solution exists. It is interesting 
to note that the condition gi(tt) = 1, which is independent of the angle 
of incidence, is the same condition on the dielectric constants that 
yields an antireflection coating at normal incidence. 

At least two layers, I ^ 2, are required to satisfy eq. (106). For 
1 = 2, e r = 1 and for a given finite 0, eq. (106) gives the unique 
solution. 

ei = [ e i + s in 2 d + (e. - sin 2 0)* cos 0]/(e* + 1), (107a) 



e 2 = eie 



(107b) 



For I > 2, the solution of eq. (106) is not unique, and we have extra 
degrees of freedom in choosing the dielectric constants of the layers. 
However, as pointed out by Baumeister, 9 the values of dielectric con- 
stants available in practice will often force a compromise of the 
effectiveness of the antireflection coatings. 

IX. DESIGN EXAMPLES 

For the purpose of demonstration and comparison of various solu- 
tions, let us design polarization-independent, multilayer, dielectric 
beam splitters with a power reflection coefficient R = 0.5 at 50 GHz, 
i.e., X = 6 mm. (The frequency is relevant only for the selection of 
the dielectric materials.) Such 3-dB beam splitters may be used in 
circular- waveguide, channel-dropping filters. 7 

All the solutions considered are symmetric, and all, with the excep- 
tion of the last solution, are of Type A with air as the surrounding 
medium; i.e., 

er= 1, 

and thus, from eq. (1) and Table II, it follows that 

k, = c ^ 1, i = 1,2, ■ ■ -,n. 

In three of the solutions, the frequency response of | r* | 2 and | r e | 2 , in 
dB, and the phase difference 

A<f> = phase (r A ) — phase (r e ) (mod 180°) 

= phase (fa) — phase (i e ) (mod 180°) 

are given. The second equality follows from losslessness and symmetry 
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Table III — Some low-loss dielectrics at millimeter wavelengths 



Dielectric 


Dielectric 
constant. 


Loss tangent 


Frequency 

range of 
measurements 


Reference No. 


Kearfott 

High-Purity 

Alumina 

Polystyrene 

Teflon 

Eceofoam PS 


9.4 

2.54 
2.1 

1 .02 to 2.0 


0.00017 

0.0012 
0.0005 

0.0004 


14 to 50 GHz 

10 to 25 GHz 
50 to 70 GHz 

10 GHz 


24 

25 

Unpublished 
report 

26 



* This is a C0 2 foamed polystyrene whose density can be adjusted to yield any 
dielectric constant in the range specified. It is not clear how the Rayleigh scattering 
by the C0 2 bubbles will affect the loss tangent at frequencies higher than 10 GHz. 



since, under these conditions, the difference between the phases of 
r and t is ±90° at all frequencies as can be deduced from eqs. (8) 
through (10). 

The dielectric materials employed in our solutions will be chosen 
from those given in Table III. Because of their low loss, these materials 
are suitable for use at millimeter wavelengths. Other suitable dielec- 
trics can be found in Refs. 24 and 25. 

9.1 Symmetric three-layer, Case 1A solution {e lf e», ei) 

With R = 0.5, eqs. (33) and (34) require that 1.64 ^ ei < 1.71 and 
e 2 ^ 15.74 for 6 = 45°; 1.46 ^ ei < 1.5 and e 2 ^ 12.34 for 6 = 60°; 
or 1.23 ^ ei < 1.26 and e 2 ^ 8.90 for = 75°. It is difficult to realize 
a solution at millimeter wavelengths for less than about 60° because 
of the high required value e 2 . However, for Fzi 75°, the value of e-> can 
be realized by Kearfott Alumina (e = 9.4). The corresponding low 
value of ei can be obtained by using Eceofoam PS, or foamed rubber 
such as Buna S Rubber whose dielectric constant can be extrapolated 
to be 1.26 at 50 GHz. 25 Thus, with ei = 1.26 and e 2 = 9.4, eq. (31) 
gives = 73.43°, eqs. (4a) and (30) give di/X = 0.428 + nn X 0.856, 
and eqs. (4a) and (31) give r/ 2 /A„ = 0.078 + m 2 X 0.172, where nn 
and m« are integers. 

The frequency response of this solution is given in Fig. 4 for mi = /»■_> 
= 0. The large required value of 6 makes this solution undesirable. 

9.2 Symmetric three-layer, Case 2A solution (e u e 2 , e,) 

With R = 0.5, eq. (46) requires that > 65.53°; and eqs. (40) and 
(42) through (44) require that ei > 4.61 and e 2 = e? > 21.23 for 
8 = 70°; ei > 2.09 and e_> = e\ > 4.38 for = 75°; and ei > 1.58 
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0.9 



0.95 1.0 

NORMALIZED FREQUENCY f/f„ 



1.05 



Fig. 4 — Frequency response of the solution in 9.1 [three-layer, Type 1A solution 
with = 73.43°, e, = 1.26 (Buna S Rubber), « 2 = 9.4 (Kearfott Alumina), rf,/\ = 0.428, 
and d 2 /\ = 0.078]. 

and e-> = e\ > 2.49 for = 78°. In the last case, the solution can be 
realized by Eccofoam PS (e = 1.59) and polystyrene (e = 2.54). Just 
as in the solution in Section 9.1, the required large value of 6 makes this 
solution undesirable. For symmetric multilayer, Type-2A solutions with 
more than three layers of two dielectrics, eq. (86) shows that d is still 
required to exceed 65.53°. 

9.3 Symmetric five-layer, Case 1A solution (d, es, &, e», d) 

For R = 0.5 and I = 2, i.e., n = five layers, eqs. (79), (80), and (81) 
can be satisfied with «i = e 3 = 2.1 (Teflon), e 2 = 9.4 (Kearfott 
Alumina), and 6 = 46.91°, which is a convenient angle of incidence. 
Further, eqs. (4a) and (67) give di/X = 0.200 + m x X 0.399 and 
d 2 /Xo = 0.084 + )« 2 X 0.168, and eqs. (4a) and (80) give d 3 /\ 
= 0.103 + >»3 X 0.399, where m x , m->, and m 3 are integers. The fre- 
quency response of this solution is given in Fig. 5 for m\ = m 2 = 
m 3 = 0. 

9.4 Symmetric seven-layer, Case 1A solution (e x , e 2 , *>, e„ e., e 2 , a) 

For R = 0.5 and I = 3, i.e., n = seven layers, eqs. (79), (80), and 
(81) can be satisfied with ei = tz = 2.54 (polystyrene), e 2 = e 4 = 9.4 
(Kearfott Alumina), and 6 = 47.52°. Further, eqs. (4a) and (67) give 
di/Xo = d 3 /\ = 0.177 + mi X 0.354 and d-,/\ = 0.084 + m, X 0.168, 
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and eqs. (4a) and (80) give d t /X = 0.026 + m 4 X 0.168, where mi, 
m 2 , and nu are integers. This solution has a narrower bandwidth than 
the solution in 9.3. 

9.5 Symmetric seven-layer, Case 1A solution obtained by cascading two 
identical, symmetric three-layer, Case 1A solutions 

(ei, e 3 , Cta 1.0, c,, e 2 , e t ) 

With R = 0.5, eq. (95) shows that R x , the power reflection coeffi- 
cient of each cascaded three-layer section, should not be less than 
0.1716; thus let us choose R z = 0.1716. In this case, eqs. (31), (33), 
and (34) can be satisfied with ei = 1.6 (Eccofoam PS), e> = 9.4 
(Kearfott Alumina), and 6 = 45.30°. Further, eqs. (4a) and (30) give 
di/K = 0.239 4- mi X 0.478 and eqs. (4a) and (32) give d 2 /K = 0.038 
+ m 2 X 0.168. 

To find the width d of the air separation, one should first find the 
phase 0x of the reflection coefficient of each cascaded three-layer 
section. By substituting the values of ei, e->, d\, di, and 6 obtained above 
in eqs. (2) through (10), one finds that <f> x = 43.59°. With R = 0.5 
and R x = 0.1716, eq. (93) gives ^ = ±90° (mod 180°), and hence, 
the electrical length of the air spacing between the two three-layer 
sections is found from eq. (92) to be <f> = 133.59° + m X 180°, where 
m is an integer. Thus, from eq. (4a) with e r = k, = 1, d/\ = 0.528 



-2.0 



-2.5 




NORMALIZED FREQUENCY f/f. 



Fig. 5 — Frequency response of the solution in 9.3 [five-layer, Type 1 A solution with 
= 46.91°, e, = « 3 = 2.1 (Teflon), « 2 = 9.4 (Kearfott Alumina), d,/\o = 0.200, 
d t /\ = 0.084, and d 3 /\ = 0.103]. 
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NORMALIZED FREQUENCY f/f n 



Fig. 6 — Frequency response of the solution in 9.6 £ three-layer, Type 2C prism solu- 
tion with = 45°, e r = 2.54 (polystyrene), e, = 1.59 (Eccofoam PS), e 2 = 1, d,/X = 0.439, 
and d 2 /\ n = 0.236]. 



+ m X 0.711. This solution has a narrower bandwidth than the solu- 
tion in 9.3. 

9.6 Symmetric three-layer, Case 2C prism solution (e,, &, 1.0, &, e,) 

Choosing = 45° and e 2 = 1.0, eqs. (16), (40), (50) and (51) give 
1.414 < ei < 2 and 2 < e r = e? < 4. Thus, the solution can be ob- 
tained with ei = 1.59 (Eccofoam PS) and e r = e\ = 2.54 (polysty- 
rene). Further, eqs. (4a) and (39) give d x /\ = 0.439 + m X 0.879, 
where m is an integer, and with R = 0.5, eqs. (4a), (4b), and (49) 
give do/\ = 0.236. The frequency response of this solution is given 
in Fig. 6. for m = and under the assumption that the outer surfaces 
of the prisms are perfectly matched. As mentioned at the end of Sec- 
tion V, the matching can be accomplished by using a quarter-wave 
layer having a dielectric constant ei = 1.59 (Eccofoam PS). 
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